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Abstract. 

We investigate the effects of baryon chemical potential /i on the shear viscosity 
coefficient r/ and the viscosity to entropy density ratio rj/s of a pion-nucleon gas 
mixture. We find that r\ is an increasing function of T and /1, while the ratio rj/s 
turns to a decreasing function in a wide region of T-/i plane. In the kinematical region 
we studied, the smallest value of y/s is about 0.3. 

1. Introduction 

Small values of the shear viscosity of a hot QCD matter inferred from the RHIC data 
have lead to a new concept, "strongly-interacting QGP (sQGP)". It is thus quite 
important to understand how such small values of the viscosity are realized in the QCD 
matter. At present, there are two conjectures about the behavior of the viscosity. From 
the analysis of the AdS/CFT correspondence, it has been conjectured that there would be 
a lower bound in the "shear viscosity coefficient to the entropy density ratio" rj/s > 1/Aix 
PQ. The lowest value n/s = 1/Ati ("the KSS bound") is satisfied by several super Yang- 
Mills theories in the large N c limit (strong coupling limit), which suggests that the 
bound could be universal. The second conjecture which is expected to universally hold 
is based on an empirical observation seen in many substances: The ratio n/s will have 
a minimum at or near the critical temperature [2] (see also [3]). More precisely, the 
ratio shows a gap at T c for the first order transition, while it has a convex shape for 
the crossover with its bottom around the (pseudo) critical temperature. Recall that the 
phase transition in QCD is most probably crossover at least for low densities. Therefore, 
what we naturally expect is that the ratio n/s in QCD will have the minimum at T ~ T c , 
and the numerical value at that point will be close to the KSS bound n/s ~ 0.1. 

These considerations motivated us to investigate the shear viscosity in QCD from 
the hadronic phase T ^ T c . Notice that we can indirectly study the properties of sQGP 
from below T c because physical quantities such as the ratio rj/s will be continuous at T c 
for the crossover transition. Moreover, inclusion of nucleon degrees of freedom enables 
us to investigate the dependence of n/s on the baryon chemical potential fi and thus to 
study the behavior of n/s in a wide region of the phase diagram. In this proceedings, 
we only give the outline of our analyses and show a few numerical results. More details 
and comparison with the results in literature are available in Ref. [I]. 
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2. Theoretical framework: relativistic quantum Boltzmann equations 

We compute the shear viscosity coefficient rj of a pion-nucleon gas mixture by solving 
relativistic Boltzmann equations which contain binary scatterings in the collision terms: 

-Lp^r(x, P ) = c^[f\n+r N \rj N } , a) 
p 

^p^f N (x,p) = C NN [f N , f N ] + C N *\J N , r] , (2) 

where f n,N (x,p) is the (iso-spin averaged) one-particle distribution of pions or nucleons, 
Ep' N = ^m^ N + p 2 and is the collision term representing binary (2 — > 2) scattering 
between particles i and j with the effects of statistics included. Also included in the 
collision terms are the scattering amplitudes, for which we adopt phenomenological 
amplitudes fitted to the experimental data of elastic scatterings in the vacuum. Fit 
was performed up to scattering energy -^s = 1.15 GeV, 2.00 GeV, and 2.04 GeV for 
the tot, ttN, and NN scatterings, respectively [I]. Note that the phenomenological 
cross sections are largely different from those of the low energy effective theories where 
p-meson and A resonances are not included. 

In order to solve eqs. (1), (2) which are nonlinear with respect to / 7r (x,p) and 
f N (x,p), we consider small deviation from thermal equilibrium at temperature T and 
baryon chemical potential [i. Namely, we linearize the equations for small deviations 
5f n,N as defined by f n,N = f^' + 5f n,N with thermal distribution f^' (Chapmann- 
Enskog method). The shear viscosity coefficient rj is given as a function of f$ ,N and 
5f n,N . Notice that rj can be decomposed into contributions from pions and nucleons: 

V = V* + V N , (3) 
where rf [r\ ) is given by pion (nucleon) distribution alone: rf = rf[fQ, Sf n } and 
rj N = r] N [f^ , 5f N }. On the other hand, we compute the entropy density in the thermal 
equilibrium, thus it depends only on /q' . 

Lastly, let us briefly discuss the range of validity of our framework. Our calculation 
has limitation in two different aspects. First, since we use the phenomenological 
amplitudes which are fitted to the data up to some finite values of scattering energy, we 
have to be careful if our results do not contain significant contributions from outside of 
the fit regions. This may be specified, for example, by a condition (s) + £ < s max 
where (s) and S = yj (s 2 ) — (s) 2 are the average scattering energy squared and its 
standard deviation, and s max is the maximum energy squared of the fit. Second, since 
our framework is based on the Boltzmann equations which are only justified for dilute 
gases, the density of particles must be small enough. This is achieved when A>cl where 
A is the mean-free path A = 1/na with a being the cross section, and d is the interaction 
range d ~ l/m n . Combining these two conditions, we find that our framework should 
give a reasonable description in a wide region of the hadronic phase on the T-fi plane: 
The boundary is given by (a quarter of) the elliptic curve connecting (T, /i) ~ (130, 0) 
and (0, 950) in unit of MeV. This is highly contrasted with the low energy effective 
theories whose range of validity is quite narrow. 
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Figure 1. Left: T dependence of 77 at u = 300, 500, 700 MeV (compared with the 
results of low energy effective theories). Middle: \x dependence of 77 at T = lOOMeV, 
and its decomposition 77 = 77" + rj N . Right: The ratio 77/s as a function of temperature 
at u = 300, 500, 700 MeV. 

3. Numerical results: \i dependence of 77 and 77/s 

Fig. [1] shows our numerical results of r] and r//s of a pion-nucleon gas mixture. The left 
panel is the T dependence of 77 at three different values of \x. In the range of temperature 
shown here, 77 is an increasing function of T. On the other hand, the result of low energy 
effective theories is a decreasing function of T, which is however not trustworthy because 
the upper limit of temperature where the low energy effective theory is valid is below 
T ~70 MeV. If one looks at the window 80 MeV < T < 130 MeV, one finds that 
rj increases with increasing /j. Mechanism of increasing 77 can be understood by the 
inspection of the middle panel where the \i dependence of the total rj as well as each 
contribution rf and r] N is plotted. Since 77 will be inversely proprtional to the cross 
section, one naively guesses that the inclusion of nucleon degrees of freedom will reduce 
the viscosity (because the effective cross section will enhance). This is indeed the case 
for rf . But in fact the contribution of nucleon viscosity itself is large, and thus the 
total 77 increases with increasing /x On the other hand, due to rapid growth of entropy 
density, the ratio 77/s turns to a decreasing function of T and \x in a wide region on 
the T-/x plane (right panel). In the kinematical region we investigated T < 180 MeV, 
H < 1 GeV, the smallest value of 77/s is about 0.3, which is realized at the edge of the 
validity region: T ~ 150 MeV and fi ~ 940 MeV. Therefore, with increasing T or \i, the 
ratio r//s becomes as small as the conjectured bound 77/s = 1/4tt ~ 0.1, but still keeps 
above the bound within the region of validity of our framework. Notice that the smallness 
of 77/s in the hadronic phase and its continuity at T ~ T c (at least for crossover at small 
fi) implies that the ratio will be small enough in the deconfined phase T ^ T c . 

4. Numerical results: Approaching phase boundary 

We do not expect we can accurately describe phase transitions within the framework 
of (standard) Boltzmann equations which are appropriate only for dilute gases. Still, 
with the help of the conjectures (in particular, the second conjecture) discussed in 
Introduction, we can extract some qualitative information about phase transitions from 
the extrapolation of our results towards critical T or fi. Below, we consider two cases: 
(i) towards higher T at zero /1, and (ii) towards higher /1 at low T. 
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Figure 2. Left: Global structure of [i dependence of rj/s at low temperature 
T = 10 MeV. Right: The valley structure of rj/s around ji — 950 MeV changes 
depending on the temperature T = 5, 10, 15 20 MeV. 

4-1. Towards higher temperatures - chiral phase transition 

The second conjecture tells us that rj/s will have a characteristic structure around the 
phase transition point: it has the minimum at the critical points. This implies that, in 
the hadronic phase, we will see the left hand side of the valley. This is consistent with the 
numerical results as shown in the right panel of Fig. 1. If this is true, then one can guess 
the position of the critical temperature from the curve of rj/s in the hadronic phase. 
For example, we can approximate the curve by a quadratic function of (T — To) with To 
being the reference temperature. Then, the critical temperature may be identified with 
the temperature where the slope of the curve is zero. We have done this for a pion gas 
(/i = 0) and for T = 140 MeV, and obtained a reasonable value T c ~ 173 MeV. 

4-2. Towards higher densities - nuclear liquid-gas transition 

At relatively high temperature T ~ 100 MeV, the ratio is a monotonically decreasing 
function of \i. However, there emerges a nontrivial structure at low temperature and 
at around normal nuclear density, as shown in the left panel of Fig. 2. According to 
the second conjecture, the valley structure implies the existence of the phase transition. 
Note that the valley locates at low T < 20 MeV and at high // ~ 950 MeV, which indeed 
coincides with the region of the nuclear liquid-gas phase transition. As temperature is 
increased, a critical line separating a nucleon gas phase and a nuclear matter (liquid) 
disappears at around T ~ 15 MeV, and above that temperature, there is no distinction 
between a gas and a liquid. This seems to be consistent with the disappearance of valley 
structure with increasing temperature as shown in the right panel of Fig. 2. The similar 
conclusion was obtained from the results of low energy effective theories [5]. 
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